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Astronomers recently discovered a non-trivial radial motion of stars within 
our galaxy and were able to reproduce the results using n-body numerics. 
I modeled a protoplanitary-like disk within a newtonian oo-body problem 
and considered the stability of local trajectories giving a qualitative expla- 
nation to the observed radial motion. As a prospect there will be a short 
estimate on the impact on spiral arms. 



Introduction 



Observational and numerical data ( Roskar et al. ( 2008h ) stress the dynamical con 



sideration of spiral arms as for stars have been observed to migrate radially through the 
spiral arms leaving significant distances between their observed positions and birthre- 
gions. It is well known that spiral arms are trajectories of local energetic minima with 
a very small tangental increase of potential energy thus making it possible to migrate 
long distances with low effort. This is why local barycenters migrate along these tra- 
jectories. The focus of this paper will now be the dynamics of objects in regions of 
these local barycenters. 

The model 

The model contains of two parts: the local gravitational impact and the non-local 
gravitational impacts. To consider the local impacts I assume the disk to be hot enough 
in order to take a region with radius R of a dense region or object (a star or planet to 
be build) so that the region will be nearly homogenous (density g, mass TO2 = 47r | fi ) 
in all three dimensions. The potential U then integrates 

AU = 47T7p (1) 



1 



Obviously this leads to 
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Demanding \Ui(r = 0)| < oo leads to A = 0. Further there shall be 
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Thus the force will be 
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Taking other forces into account the force will slightly change, hence will be pertur- 
bated to 

Fi = -a 2 r , a 2 >0 (6) 



Integrating all the non-local impacts you will get a barycentral force 

F e = — i , ai > 
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With ai :— m being the mass of the dense object, Newton's second law becomes 

(8) 
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Whereas (x; ip) T is the current position of the local barycenter with respect to the 
barycenter of the non-local masses and (x; y) T being the position of the dense object 
with respect to the barycenter of the non-local masses being in the center. 

I do not consider any motion normal to the plane of the disk, since it is obvious that 
it will be just some oscillation. The model also only applies to the middle regions of 
the disk, since there are non-trivial ^-components to the force inside the bulge and in 
the outer regions the considered region will not be as homogenous anymore. 

Stability of trajectories 

Now the trajectories can be characterized by 
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with X :— (x, v x ,y, v y ) T and x 2 + y 2 ^> R 2 . Obviously is / in X lipschitz continuous 
and thus complies with the standard prerequisites. A possible Lyapunovfunction will 
have to integrate 



0>F:=(W|/)= v x d x V + 
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Using the identities x 
X' := (x>,v' x ,y',v' y ) T 
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integrates [10] in the open region 

n := <! (x,v x , y,v y ) T : \x-x\ < 
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with V(x, v x , ip 7 V^) — = V(x, v x , ip, V^) and coordinates respectively chosen to 
comply x^y^VxjVyjXjip > 0, which is perfectly possible since x,y,x,i> > can be 
chosen as there is x 2 + y 2 3> R 2 > (x — x) 2 + (y — V") 2 an d v x ,v y > define the 
surface normal. The implicit function theorem ensures that there will always be a 
time-interval to locally solve these conditions. Should the Jacobian be singular then 
simply turn the coordinates a little (a; 2 + y 2 ^> R 2 > (x — x) 2 + {y — ip) 2 )- ^ herein is 
the largest possible local region to be observed. 
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Obviously 
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thus in every region of (x, v x , tp, ity) exists a point Z with V(Z) < 0. Hence the local 
region is instable and as a consequence any dense object within the region with leave 
this in a finite period of time. 



On migration 

Since it is known that the local barycenters travel along the spiral arms any dense 
object in that regions will move with them and this way migrate radially through 
the disk. On the other hand each object will leave that region after a finite time 
and hence travel no longer with that local barycenter. Now the object is no more in a 
local energetic minimum and hence will move back entering another radially migrating 
region just to soon leave this again and enter right another. Multiple objects of the 
same birthregion will do so indepenently. Thus the observed galaxies should show a 
good mixture of young and old stars and stars with different chemical compositions 
since there are stars from basically every birthregion close another. These are exactly 
the observed properties in the middle parts of galaxies. 

It might be interesting to estimate the density of stars in spiral arms for long term 
effects. In order to migrate outwards for an object it needs energy which has to be 
taken from other objects for the sum of energies to stay conserved. Since it seems 
more likely to transfer low amounts of energy I will approximate the energy-transfer- 
distribution as gaussian with the expected transfer amount [Ie = 0. Let then T be 
the average time between entering two consecutive regions and let A be the average 
radial movement, then the energy-transfer-distribution ipE with standard deviation oe 
becomes 

1 / A.E 2 \ 

In order to estimate oe I considered the energy E g — ai 2 " 3 of a mass m 3 at radius r. 
As the radius changes by A the transfered energy is given by 

If |ai|m 3 |ai|m 3 \ _ |ai|m 3 / 2r \ 
^ E - 2\2{r + A) + 2{r - K)) - \t^~^) (15) 

Now I average this over the entire applicable region with inner radius i?j and outter 
radius Ro 

1 f Ro , \ai\m 3 , fR 2 n -A 2 \ , . 

^=R^R-J Rl ^ Edr= A(R^RV) ln {-Rj^) (16) 



4 



Now it is possible to approximate oe by 

1 /"Mae 

- = 2 / <p E dAE (17) 
z Jo 

The energy-density thus will be homogenous. Considering our own galaxy with Ri = 
8000Zy < r < 50000/y = i? one finds E(8000ly) w 6ail0~ 21 ^m 3 as well as 

£(50000Z|/) « ail0- 21 ^m 3 and thus A£ = 5Qil0^ 21 ^m 3 on a distance of 4200(% 
which is nearly homogenous. 

Conclusion 

From (fl5|) it is easy to see that the energy-transfer enables much larger distances 
of migration the farther away an object is from the center of the disk. Hence one will 
more likely find a higher mixture of objects in the outter regions. On the other hand 
this gives a mechanism preventing too many masses to fall into the interior regions 
and thus from possibly falling into black holes. Further keeps the disk from accreting 
too fast making it possible for planets and stars to form. 

The local instability hence has an impact on the largescale stability of a protoplan- 
itary disks and spiral arms. It does also explain the observed mixture of stars. 
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